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Introduction 


The  best  single  book  on  the  subject  under  review  'ferf-  contains  some 

»»  no,  urU.ft'i, 

350  references  and  is^*ew  over  5 years  old.  The  method  adopted  in  this 
report  is  to  start  with  a very  brief  outline  of  the  theory  in  that  monograph  by 
Beckmann  and  Spizzichino,  specialize  to  the  case  of  scatter  in  the 

1 r 

specular  direction  and  bring  the  theory  up  to  date  (Section  |) . Section(IIj 
provides  simplified  prediction  equations  and  graphs  based  on  the  theory  of 
Section  I, 

SectionvIII,ds  devoted  to  a summary  of  the  major  known  inadequacies 
of  the  present  prediction  scheme  and  brief  statements  of  problems  that  must 
be  solved  before  we-can^valuate-and-improve  the  accuracy  of  the  current 
predictions  v ^ ° <■*>/  .*  iv  * I <*  > 1 ■>  « i n,  y f<w  c <1 0 


AfiST,-  MC  (, 


Theory 


A.  The  Fully  Developed  Sea 

We  assume  that  the  sea  surface  from  which  the  sound  is  scattered 
is  homogeneous,  and  that  the  temporal  variation  of  sea  heights  at  a point 
or  the  spatial  variation  at  any  instant  of  time  [2]  are  Gaussian  distributed 
with  mean  value  zero, 


m 


and  with  variance , a , 


<C2>  - a2 


i 1 

V ’*-• 
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The  surface  height  correlation  function  is  C 3 ] 

c , Z.  = <C  (0.0,0 ) Ctu.VjTL) 

a2  o 2 

and  the  Fourier  transform  of  Z yields  the  three-dimensional  spectrum  of 
the  sea 


X(xx,x  n).(2ir)'3^  Z(u,v,t 

y ‘TJU,V 


) exp  ti(x  u + X v - fi  t)  ] dudvdT 
x y 


where  u 
v 

r 


x ,x 
x y 

0 


= displacement  in  x direction  between  two  surface  points 
= displacement  in  y direction  between  two  surface  points 

* displacement  in  time 

* components  of  surface  wave  propagation  constant 
= angular  frequency  of  surface  wave 


We  shall  assume  the  gravity  wave  relation  between  angular  frequency, 
propagation  constant  and  acceleration  of  gravity,  g: 

2 

0 = xg 

This  permits  the  spectrum  of  the  sea  to  be  specified  in  two  dimensions 
either  as 


X = * (x  , x ) 6 (Vgx  - Q) 
x y 


or 


X = * (x  ,a)  6 (Vgx  - 0) 
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Jjr 


-m 


where  the  Dirac  delta  function,  ft  (vgx  - 0) , simply  selects  the  wave  that 
fulfills  the  gravity  wave  relation. 

The  more  commonly  measured  frequency  spectrum  $»(0)  is  now 
obtained  by  using  the  energy  density  coordinate  transformation 

'I'xdx  = <t>  (0)  d ft  . 

The  frequency  spectrum  is  more  simply  expressed  as  the  integral  of  the 
polar  form  of  the  x spectrum  over  all  azimuth  angles 


$(fi) 


2 a 


r.3  2n 

Y \ ^(x,  a)  da 


(1) 


The  semi-empirical  form  of  the  frequency  spectrum  for  the  fully- 
developed  sea  that  is  most  widely  accepted  today  (Pierson- Moskowitz)  [4] 
is 


*(0)  = 


m 


Q" 


exp 


MW] 


(2) 


where 


a 

J8 

O 

c 

W 


19.5 


= 8.1  x 10’3 
*=  0.74 

S g/^19.5 

=s  Wind  speed  at  19. 5M  above  sea  surface 

=*  acceleration  of  gravity  in  units  consistent 
with  those  of  W . 
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The  three  parameters  of  the  sea  that  are  most  significant  in  specifying 


the  sound  scatter  are  the  mean  square  surface  height,  o , the  mean  square 
2 

surface  slope,  £ , and  the  surface  correlation  function,  C. 


We  obtain  the  mean  square  height  by  integrating  the  frequency 
spectrum,  Eq.  (2): 


cu 

a2  = J 0(0)  dfl 


For  the  mean  square  slope,  we  will  use  the  relation  determined  from 
optical  measurements  of  the  sea  [5] 


5.12  x l(f  5W.,  + 0.003  + 0.004 
41  — 


where  is  measured  in  cm/sec  at  a height  of  41  feet  over  the  surface. 


Unfortunately  the  measured  speed  is  sensitive  to  the  height  of  the  anemometer. 


Finally,  for  simplicity,  and  because  of  lack  of  sufficient  two- 
dimensional  sea  data,  we  will  assume  that  the  spatial  correlation  function 


is  of  Gaussian  form  and  isotropic: 


C (1)  = 


where 


1 V 2 a.  2 

1 = *U  + V 


and  L = surface  correlation  length. 


The  assumption  of  an  isotropic  Gaussian  function  makes  it  possible  to 
calculate  the  correlation  length  from  a knowledge  of  rms  slope  E and  rms 
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height  a by  the  relation  that  is  valid  for  seas  of  lesser  roughness  Cl] : 


E 


(6) 


B.  Sound  Scattering  Theory  (Isakovitch-Eckart-Beckmann) 

The  theoretical  solution  to  the  sea  surface  sound  scattering  problem 
usually  starts  with  the  Helmholtz  Integral: 


PQ 


4"  I r2 


dS 


where 

S 

k 

X 

r2 

n 


the  scattered  acoustic  pressure  at  an  interior  field  point,  Q, 
of  an  enclosed  volume. 

integrating  surface  selected  so  that  acoustic  field  is  zero 
everywhere  except  on  the  insonified  water  surface. 

2ir/\  = acoustic  propagation  constant 
acoustic  wave  length  in  water 

distance  from  surface  scattering  region  to  field  position, 
normal  to  surface  scattering  point. 


Fig.  1 Scattering  Geometry 


-5- 


After  applying  the  boundary  conditions  for  acoustic  pressure  and 
acoustic  particle  velocity  at  the  sea  surface,  and  assuming  that  there  is  no 
shadowing  or  secondary  scattering  even  for  near-grazing  incidence,  we 
obtain  the  integral  for  the  scattered  pressure  from  a given  surface  of 
scattering  area  AA  (now  called  Ap2  H Pq)  in  terms  of  the  geometry: 


Ap2  = Ap2.  h L*lkl'v  '~dS 

plane  AA 


where  Ap2  = scattered  pressure  for  a plane  surface  area  AA. 

plane 


A A 


*1 


~2 


* ikp.  AA  cos  0]  e^^r2 

- insonified  surface  area  over  which  p^  and  0^  are 
constant. 

= incident  sound  pressure  at  sea  surface. 

= incident  sound  angular  frequency 


1 + cos  0^^  cos  ©2  - sin  ©1  sin  ©2  cos  0^ 
cos  0^  (cos  01  + cos  02) 


= angle  of  incidence  (measured  with  normal) 

= angle  of  scatter  (measured  with  normal) 

= azimuthal  angle  of  scatter  plane  with  respect  to 
incident  plane. 


= k (ij  sin  ©1  - lg  cos  6^ 

= k (j^  sin  02  cos  03  + sin  ©2  sin  ©3  + cos  d^) 
il , i_2 'ig  are  the  unit  normals  in  the  x,  y and  z (depth)  directions, 


respectively. 


es»'i 


^ ' '~^- "i^^SsfcSF 


The  mean  square  scattered  pressure,  which  in  general  includes  both 
coherent  and  incoherent  sound  intensity,  is  formulated: 


<Ap2Ap2*>  . Ipl  -E-J  $$$$  (eiKzCC<x.y)-C(x*,y')l)  0i  [K^u  + Kyv  3 

*V  I i 

x y x'y 


(dxdx’dydy1) 


where 


S.  ’ h'h 


U = X*  - X 


= y - y 


h Ap2* 

plane  plane 


= complex  conjugate 


= ensemble  average 


The  ensemble  average,  < ) , can  be  easily  evaluated  for  our  surface 


with  its  two-dimensional  Gaussian  distribution: 


<eiKz(C-C'))  a e-Cko(cos  61  + cos  02)]  [l-C] 


Define  the  ACOUSTICAL  ROUGHNESS  OF  THE  SEA  SURFACE  a R 


R 2 Ckcr(cos  8^  + cos  82)] 


so  that  the  mean  scattered  intensity  can  be  written 


<ap2ap2«>  = e~*^~c)  $$n el(K*u + r'yV>  d5K,x'dydy' 


x y x'  y' 
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We  are  interested  in  the  relative  scattered  intensity,  (compared  with  the 
plane  surface  "scattered"  intensity)  specialized  to  the  requirements  of 
this  report,  specular  scatter,  in  which  0^  = 9,  , 0^  = 0 and  F = 1. 


4*-. 


<APaP*)  : 

<«>*)  « — e * ~A 

plane 


(e-R  E 
V m=l 


The  two  terms  of  Eq.  (8)  represent  coherent  scattered  sound  (the 

~R 

e term)  and  incoherent  scattering  (the  summation  term) . A graph  of  the 
term  in  parentheses  (Fig.  2),  defined  as 


2 e 


shows  that 


S (R)  ■»  R for  R « 1 


S (R)  * for  R » 1 


Since  it  can  be  assumed  that  L < A A,  we  can  conclude  that,  in  the  limits, 
the  specular  scatter  for  a Gaussian  sea  will  be: 


For  Small  Acoustical  Roughness,  R <<  1:  (pp*>  - e 


For  Large  Acoustical  Roughness , R » 1:  ^PP*)~^TT- — r” 

KAA  2k  E AAcos  0j 


It  is  important  to  observe  that  the  value  of  db  loss  per  bounce  is  an 
inadequate,  and  in  fact  misleading,  concept  for  the  rouvh  surface  case. 
In  that  limit  (R»  1)  the  specific  value  of  the  scattering  area  becomes 


significant  and  the  asymptotic  value  of  <pp*>  will  change  from  one  bounce 
«.o  the  next  as  the  sound  wave  continues  to  diverge  and  scatter. 

In  addition  it  should  be  noted  that  for  an  isotropic  Gaussian  sea 

the  specular  scattering  loss  in  the  low  acoustical  roughness  limit  is 

dependent  on  the  rms  surface  height  and  is  completely  determined  by  the 
2 2 2 

value  of  R = 4k  a cos  6^  . On  the  other  hand,  in  the  large  acoustical 
roughness  limit,  the  scattering  loss  factor  will  vary  inversely  as  the  mean 
square  surface  slope  (as  well  as  the  acoustic  propagation  constant,  the 
angle  of  incidence  and  the  scattering  area) . 


C.  Some  Other  Sound  Scattering  Theories 

The  theory  outlines  in  Section  IA  had  its  origins  in  work  by 
Isakovitch  (1952)  and  Eckart  (1953) . A different  approach,  based  on  a 
generalization  of  Rayleigh's  solution  for  a corrugated  surface  was 
presented  by  Marsh  et  al. , [?]  and  criticised  as  incorrect  by  Uretsky 
and  others.  Rather  than  review  the  theory  and  the  objections  of  its  critics, 
the  appendix  of  this  report  will  point  out  the  inadequacy  of  predictions  from 
Marsh's  theory.  It  will  be  seen  that  Marsh’s  theory  makes  no  provision 
for  variation  of  angle  of  incidence,  that  it  is  consonant  with  the  pre- 
dictions of  the  development  of  Section  B only  for  losses  of  zero  db  to 
approximately  3 db  after  which  it  diverges  wildly  from  the  conclusions  of 
others  and  from  experiment. 

Uretsky  [8]  has  presented  an  exact  solution  to  the  related  problem 
of  scattering  from  a pressure-release,  sinusoidal  corrugated  surface.  This 
work  has  been  used  in  a laboratory  model  study  by  Barnard  et  al. , [9]  and 
has  been  found  to  give  an  accurate  description  of  this  type  of  scattering. 
The  applied  interest  here  lies  in  the  expectation  that  scattering  from  swell 
may  be  similar  to  scattering  from  a corrugated  surface. 
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II 


Prediction  Procedures 


A.  The  Relative  Contributions  of  Coherent  and  Incoherent  Components 
The  conclusion  of  Section  IB 

<PP*>  = e"R  + S(R)  (8) 

? 2 

where  x = :r’  sin 

> <• 

A 

to  = grazing  angle  = 90°  - 
_ -m 

S(R)  = e"R  S — r (9) 

, mmi 
m®l 

L surface  height  correlation  length 

constitutes  the  best  current  means  for  estimating  the  relative  mean  square 
scattered  pressure  (compared  to  a smooth  surface)  for  specular  scatter. 

The  known  sources  in  precision  and  error,  due  -o  the  assumptions  that 

lead  to  Eq.  (8) , are  discussed  in  Section  III.  ! 

*- 

In  order  to  proceed  to  a simplified  prediction  procedure,  we  present  1 

in  Fig.  3 a "typical"  graph  of  Eq.  (8)  plotted  with  R as  independent 

| 

variable,  :based  on  the  assumptions  that  the  increase  in  R occurs  because 

2 « 

of  increasing  frequency,  and  that  there  is  no  change  in  k A A (piston-like 

i 

i 

source),  81  , or  sea  condition  (cr,  L,  2).  For  small  acoustical  roughness 

* i 

the  scattering  is  coherent  and  decreases  exponentially  as  R increases.  j 

The  incoherent  scattering  term  is  initially  very  small,  it  increases  with  R . 

and  approaches  an  asymptotic  value  determined  by  the  rms  slope  Z . ; 

1 

The  value  of  this  asymptote  depends,  importantly,  on  the  inson-  j 

| 

ified  area,  A A.  For  the  example  plotted  we  have  assumed  an  insonified  i 


11 


2 2 2 4 

area  of  sizo  in  k A A ■ 4ff  AAA  * 10  . A larger  value  of  tiro  area  will, 

of  courso,  docroaso  the  magnitude  of  tho  incoherent  component  so  that  the 

**■ 

cohoront  bohavior  e will  continue  to  dominate  to  larger  values  of  R , 

B.  Decibels  per  Bounce  for  the  Coheront  Component  Alone 

If  we  assume  that  tho  uniformly  insonified  surface  area  is  very 
2 6 

large  (k  AA  »10  ) thon,  almost  regardless  of  the  possible  sea  slope  and 
grazing  anglo,  tho  simple  cohoront  scattering  relation  <pp*>  » e~*  will 
give  the  whole  picture.  Fig,  4 shows  this  behavior,  now  expressed  in 
db/bounco,  calculated  from 


db/bounco 


10  log1Q  e 


as  a funotion  of  grazing  angle,  wind  speed,  and  rms  height  of  the  sea  for 
tho  frequencies  50,  100,  800,  1200  and  3500  Hz  and  the  grazing  angles 
» 5° , 10°  and  20° , basod  on  tiro  simplifying  assumptions  that  have  been 
made.  The  Piorson-Moskowitz  [4]  spectra  for  a fully-developed  sea  have 
been  used  to  determine  a from  the  wind  speed.  The  flaws  in  these 
assumptions  and  tho  task  wo  face  before  we  can  improve  the  prediction 
accuracy  will  be  discussed  briefly  in  the  next  section.  The  limited 
experimental  support  for  those  theories  is  summarized  in  the  Appendix. 
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Fig.  4 DB/Bounce  for  Coherent  Scattering  from  a Pierson-Moskowitz  Sea. 
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Future  Prediction  Procedures.  Unsolved  Problems 


A.  Surface  Shadowing 

The  theory  of  Section  IB  assumes  that  there  is  no  surface  shadowing 
or  secondary  scatter  from  one  wave  to  another.  Preliminary  attempts  to 
correct  for  near-grazing  shadowing  C 1 1 3 have  been  based  on  a ray  acoustics 
approach.  It  is  clear  that  diffraction  into  the  shadow  region  must  be 
considered  and  that  the  large  corrections  that  Wagner  obtains  will  be 
decreased  for  the  cases  of  Navy  interest.  The  fact  that  the  correction 
should  be  made  ’•■hrn  crazing  angle  is  approximately  equal  to  or  less  than 
rms  slope  (e.g.  <p  < 6'  £t  W = 4 knots,  <p  < 13°  for  a 20-knot  wind) 
makes  it  impo'-  art  that  this  correction  be  evaluated  as  soon  as  possible. 

It  could  not  be  done  within  the  limits  of  the  present  report. 

B.  Near  Surface  Bubbles 

As  sound  incidence  approaches  grazing,  absorption,  dispersion 
and  scattering  by  near- surface  bubbles  can  become  so  important  that  the 
sound  effectively  either  does  not  reach  the  surface  or  it  changes  angle  of 
incidence  and  amplitude  in  approaching  the  surface.  The  number,  size, 
and  depth  of  entrainment  of  surface-generated  bubbles  increases  with 
increasing  wind  speed.  The  scant  literature  on  the  subject  includes  a 
laboratory  study  of  bubble  generation,  one  or  two  in- situ  measurements 
and  a few  theoretical  studies  (e.g.  ref.  [12]).  There  was  inadequate  time 
during  the  present  study  to  evaluate  the  effect. 

C.  Non-Fully-Developed  Seas 

Typically  a sea  will  show  different  rms  slopes  in  the  upwind-downwind 
direction  compared  with  the  cross-wind  direction.  This  will  affect  the 
incoherent  component  (through  E)  which  we  have  assumed  to  be  negligible 
for  the  purpose  of  simplifying  Section  IIB. 

-15- 
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If  we  are  talking  about  an  attenuated  sea , propagated  from  distant 
storms  ("swell") , the  spectrum  will  no  longer  be  the  one  described  by 
Pierson  and  Moskowitz.  The  high  freque  ..cies  will  have  dropped  out  much 
more  than  the  lows  and  the  spectrum  will  be  much  narrower.  Presumably, 
if  the  swell  spectrum  is  sufficiently  narrow,  the  behavior  will  be  similar  to 
scattering  from  a corrugation,  whereas  if  the  frequency  spread  of  the  swell 
spectrum  is  great  enough  the  Gaussian  sea  approximation  would  be  more 
accurate.  The  literature  of  scattering  theory  is  completely  silent  on  this 
subject.  We  still  do  not  know  when  to  solve  the  problem  of  scattering  from 
swell  by  using  the  solution  for  a sinusoidal  corrugation,  or  that  for  a 
Gaussian  sea,  or  some  compromise  between  these  two  extremes.  The 
importance  of  this  decision  is  clear  if  it  is  realized  that  scattering  from 
a corrugation  will  vary  from  zero  to  a maximum  as  the  angle  of  incidence 
goes  from  (a)  parallel  to  an  element  of  the  corrugation  to  (b)  in  the  plane  of 
the  sinusoid.  On  the  other  hand,  for  an  isotropic  Gaussian  sea  we  have 
seen  that  the  coherent  term  does  not  depend  on  the  plane  of  incidence. 

If  the  sea  is  not  fully-developed,  the  relation  that  yields  rms 
height,  a , from  wind  speed  will  depend  also  on  the  duration  of  the  wind 
flow  at  that  speed  and  on  whether  growth  occurred  from  a white  noise 
background  or  from  a flat  calm  ocean.  It  is  now  possible  to  estimate 
this  effect  [13]  . 
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APPENDIX 

Comparison  with  Experiment 

Verification  of  the  accuracy  of  scattering  predictions  from  Eq.  (8) 
have  been  foiled  at  sea  by  unknown  or  poorly  defined  measurement  of  the 
oceanographic  parameters.  Specifically,  use  of  Eq.  (8)  requires  knowledge 
of  the  rms  height,  cr , the  rms  slope,  £ , and  the  area  which  is  uniformly 
insonified.  If  the  acoustical  roughness,  R,  <<  1,  the  only  oceanographic 
parameter  required  for  the  prediction  of  db  loss/bounce,  is  the  rms  height. 


Looking  at  the  history  of  ocean  surface  scattering  experiments  we 


find: 


(A)  Marsh  [10]  shows  data  from  a large  number  of  unpublished  ocean 
experiments  in  order  to  compare  with  his  theory 


db/bounce  * - 10  log  [l  - 0.0234b'3'  ] 


b * (F,kcps)  (Wave  Height,  H,  in  Ft.) 


Marsh  assumes  that  there  is  no  explicit  dependence  on  angle  of  incidence. 
However,  if  a Gaussian  distribution  is  assumed,  we  can  take  H = 6or(to 
include  99%  of  the  values  of  height) , and  convert  Marsh's  formula  to 


bounce 


= - 10 


iog10  [i-t.omf^b-)372] 


Unfortunately,  no  angles  of  incidence  or  sea  conditions  are  given  for  Marsh's 
data  which  fits  his  theory  reasonably  well  as  long  as  db/bounce  < 3.0. 

Marsh's  data  for  R < 1 also  fits  quite  well  to  Eq,  (8)  if  we  make 
the  reasonable  assumption  Qj  = 80°  (10°  grazing). 


lik. 


The  theory  developed  by  Marsh  has  been  rejected  in  this  report 
because:  it  is  valid  only  for  R < 1;  it  contains  no  explicit  provision  for 
introduction  of  dependence  on  angle  of  sound  incidence. 


(B)  Project  AMOS  [14]  includes  an  estimate  of  db  loss/bounce  giyen  in 
terms  of  depth  of  the  isothermal  layer  which  can  be  restated  as 


db/bounce  = 2.3^f*  for  sea  state  < 3 

db/bounce  « 4.6VT  for  sea  state  > 3 

where  f = sound  frequency  in  kHz. 

These  formulas  are  too  crude  and  too  divergent  from  recent  studies  to  be 
considered  further.  The  data  may  have  been  affected  by  bubbles. 

(C)  Proud,  Beyer  and  Tamarkin  [15]  performed  a laboratory  study,  using  a 
one-dimensional  cork  model  of  a Gaussian  Sea,  and  showed  that  the  theory, 
such  as  presented  in  Section  IB,  is  essentially  correct,  and  that  forward 
scattering  depends  only  on  acoustical  roughness 


16ff2g2cos2&i 


slope  S and  area  insonified  for  angles  of  incidence  from  0°  to  60°,  for 
0.3  < R<  2. 

(D)  Medwin  [6]  used  a wind-driven,  scaled,  model- sea  surface  of  a large 
tank  of  water  to  study  mean  value  and  probability  distribution  of  normally 
incident,  forward- scattered  acoustic  pressure  over  the  range  0.25  < R < 25. 
The  oceanographic  parameters,  rms  height,  height  distribution,  rms  slope 
and  its  distribution  in  the  up-down  wind  direction  and  cross-wind  direction 
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were  measured  directly  in  order  to  check  the  predictions  of  theory.  Accurate 
predictions  were  obtained  when  Eq.  (8)  was  used. 


(E)  Scheible  et  al. , [16]  have  used  the  same  laboratory  set-up  as  in  (D)  to 
study  near-grazing  forward  scatter  in  the  case  II  < 0.3  where  only  the 
coherent  scattering  term  should  be  important.  Their  work  suggests  that 
near-grazing  scattering  lies  between  the  value  given  by  the  coherent  term 
of  Eq.  (8)  and  the  ray-corrected  expression  due  to  Wagner  [ll]. 


^ \ M 

• V-Sj 

i r J 


j 

p f j 
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